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Abstract 

The Curie- Weiss model is an exactly soluble model of ferromagnetism that allows one to study in 
detail the behavior of thermodynamic quantities in the neighborhood of the critical temperature. 
Because of its simplicity and because of the correctness of at least of some of its predictions, 
the Curie- Weiss model occupies an important place in the statistical mechanics literature and its 
application to information theory. It is frequently presented as an introduction to the Ising model or 
to spin glass models, and usually only general features of the Curie- Weiss model are presented. We 
discuss here properties of this model in a rather detailed way. We present the exact, approximate 
and numerical results for this particular model. 
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I. INTRODUCTION 



A phase transition is the transformation of a thermodynamic system from one phase or 
state to another. During a phase transition of a given medium certain properties change, 
often discontinuously, as a result of some external conditions, such as temperature, pressure 
or the magnetic field. For example, a liquid may become gas upon heating to the boiling 
point, resulting in an abrupt change in volume. 

In the case of a ferromagnet one should predict the dependence of such quantities as the 
free and internal energy, heat capacity and magnetic susceptibility on temperature T and 
magnetic induction B. Generally, this is very complicated task. Mar stressed the distinction 
between the direct approach to the problem of phase transitions and the approach exploiting 
symmetries of the problem. Here we shall illustrate the former approach. This means 
calculations of physical properties of interest in terms of parameters given in a particular 
model, i.e. solving a model. The calculations may be done analytically or numerically; 
exactly or approximately. 

The approach exploiting symmetries does not attempt to solve the model. From various 
symmetry properties one deduces some important physical properties. 

We shall focus on a particular model - the Curie- Weiss model of a magnet. This model 
is also called the infinite range Ising model. The Curie Weiss model is an exactly soluble 
model of ferromagnet ism that allows one to study in detail the behavior of thermodynamic 
quantities in the neighborhood of the critical temperature. Since not all predictions of this 
model agree with experiment, other models must be considered. However, because of its 
simplicity and because of the correctness of at least some of its predictions, the classical 
Curie- Weiss model occupies a central place in the statistical mechanics literature. It is 
frequently presented as an introduction to the Ising model or to spin glass models^ - —. In 
these references only general features of the Curie- Weiss model are discussed. This is why 
we shall discuss here properties of this model in a rather detailed way. We present the exact, 
approximate and numerical results for this particular model. 
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II. THE CURIE- WEISS MODEL 



One of simplest classical systems exhibiting phase transition - the Curie- Weiss model has 
been introduced by Pierre Curie and then by Pierre Weiss in their development of simplified 
theory of ferromagnetism^. We shall call it the C-W magnet. Let us call the set of integers 
from 1 to iV a lattice, and its element % a site. We assign a variable Si (the Ising spin) to 
each site. The Ising spin is characterized by the binary value: +1 if microscopic magnetic 
moment is pointing up or —1 if it is pointed down. Particles with the Ising spins interact 
via Hamiltonian 

Hint = ~ Yl s * s r W 

l<i<j<N 

The constant J is positive. The interaction energy of all pairs of spins of the Curie- Weiss 
magnet is the same and their interaction depends on N. The normalization by 1/N makes 
Hn a quantity of the order N. The underlying assumption of an infinite-range interaction is 
clearly unphysical. The Hamiltonian ([I]) does not depend on dimension of the space which 
Curie- Weiss magnet is occupying. 

The magnetic moment of a particle is proportional to the spin /Zj = /iSj, where fi is 
the magnetic moment. In an applied magnetic field with the magnetic induction vector B 
particles with magnetic moments being parallel or antiparallel to B, acquire the energy 

N 

H f = -iiB Si ■ (2) 

i=i 

The complete Hamiltonian consists of two terms 

j N 
H = ~— Yl s iSj -fiBj2 s i- ( 3 ) 

l<i<j<N i=l 

The Hamiltonian ([3]) does not change if we reverse signs of all spins Sj — > — Sj (i = 
1,2,..., N) and the direction of the induction vector B — > — B 

H (si, . . . , s n , B) = H . . . , -s n , -B) . (4) 

Denote a particular configuration (si, S2, ■ ■ ■ , sjv) by {s}. To each configuration {s} there 
corresponds an energy E ({s}) = H ({s}). 

The nature of phase transitions of magnetic systems is well understood. At temperature 
OK magnetic systems, in particular the M-W magnet, are in a lowest energy state with all 
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spins being parallel. Thus, their magnetization M is finite and our magnet is ferromagnetic. 
As temperature is increased from zero the thermal noise randomizes spins. A fraction of 
them become antiparallel. This disorder grows with raising temperature and a diminishing 
fraction of them points at the initial direction. At temperature T c - the critical temperature, 
and beyond, magnetization vanishes and the material becomes paramagnetic. For T above 
T c , there must be macroscopically large regions in which a net fraction of spins are aligned 
up. However, their magnetization mutually compensates - they cannot make a finite fraction 
of all regions agree. For T just below T c the compensation is not complete and the small, 
but finite, fraction points in the same direction. 

For left-hand vicinity of T c thermodynamic functions depend on the dimensionless pa- 
rameter t = (T — T c ) /T c and consists of terms regular in t and singular in it. The singular 
terms depend on powers of These powers are called the critical indexes (or critical 

exponents) and are defined for B = and t — >■ 0. The specific heat c per one particle and 
the magnetic susceptibility xt for the paramagnetic phase are characterized respectively by 
critical indexes a and 7, whereas in the ferromagnetic phase they are characterized by a' 
and 7'. For example in the paramagnetic phase 

c ~ \t\~ a , Xt ~ |t|~ 7 . ( 5 ) 

We shall calculate these sets of critical indexes, as well as the dependence of the internal 
energy and entropy on \t\, for ferromagnetic and paramagnetic phases of the Curie- Weiss 
magnet. 

In the presence of the magnetic field the dependence of magnetization on the magnetic 
field at \t\ — is characterized by the critical index 5 

m ~ B 1/s . (6) 



In the following, in place of temperature we will use 9 = ksT, fc B being the Boltzmann 
constant. 
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III. CALCULATION OF FREE ENERGY 



Note that due to the particular structure of the Hamiltonian of C-W magnet, the Hamil- 
tonian (CO can be written as 



H, 



J 



N 



hit 



2N\^ Sl ) + 2" 

. i=l 



(7) 



This form of Hamiltonian can be considered as a defining feature of mean field models. 
The partition function is defined as usual as 11 



Z N = ^ e" 

W 



E({s})/9 



(8) 



The summation is performed over all 2 N configurations {s}. Note that a different method 
calculating spin configurations can be used^ - -. 

Introduce two dimensionless quantities K = J/8 and h = \iBjQ. Then, the partition 
function (181) can be written as 



Z N (9) = 22 exp 



e 2 exp 



N 



N 



.4=1 

N 



2iV ^ S; 



i=i 



i=l 



(9) 



Thermodynamic functions can be obtained via the Helmholtz free energy F N (0, By± 



F N (9, B) = -9 In Z N (9, B) . 



(10) 



Free energy is minimized when a system reaches equilibrium^^. The derivative of F^ (9,B) 
with respect to B gives 



dF N (9,B) 
dB 



N 



N 



(11) 



1=1 



1=1 



where (sj) means the mean value of Sj calculated with the canonical distribution function^. 
The minus derivative (llip defines magnetization M of C-W magnet, thus 



iV 



M = /i^( Sj ). 



(12) 



i=i 



The second derivative of Fn (0, B) with respect to B is positive^ 



Since [d 2 F (9, B) /dB 2 ] e defines the magnetic susceptibility xt (cf. Eq. we conclude 

that the magnetic susceptibility is positive. Therefore, the change of induction 5B increases 
the magnet energy by 6B-X-SB > 0. The CW magnet fulfills the condition of thermodynamic 
(macroscopic) condition of stability^. 

A simple evaluation of the partition function fl9]) is precluded only by square of magne- 
tization in the exponential. One can get rid of this square using the Gaussian linearization 
of the form 



+oo 

„2 



e 



TBI*"*- (14) 



N 



In the present case a = y/K/ (2N) s 4 . Now, the partition function factors with respect 

i=i 

to individual summations over the state s,- 



_i£ +1 +1 +1 

v e E ••• E 



2tt 

Sl=— lS2=— 1 Sjv = -1 



s N £+hs N 



(15) 



N „—K 



2 n e 



die 



cosh (^v/^/iV^ + Zi 



N 



— oo 



Performing the change of variable a/ K/N£ = Ky we get 

1/2 + °° 

Z N = 27V (^) e " f / d V^K,y (V) , (16) 

— oo 

where 

®K,h (y) = e- Xy2 / 2 cosh (Jfy + /i) . (17) 

In addition to y the function $>K,h (y) depends on two dimensionless parameters K and h, 
i.e. on 8 and B. 

Free energy per one particle is proportional to (InZ/v) /N. Since we are interested in 
analyzing the system in the large size limit 

lim — = lim In ( Z\l N \ — In ( lim zU N 

JV-s-oo N JV-s-oo V / \N^oo JV 



For Zn (|T5|) we obtain 

_f(e i B) = 1 l n ( e -f JrN/2-k) + In 2+ 

l/N (lg) 



+oo 



+ ln Km { / (y)] 

iV— >oo 



N 



In order to obtain the explicit form of the function / (9, B) we use the Laplace theorem^. 
The result is 

- ^ -In max $ AVl (y) + ln 2 . (19) 

(7 — oo<j/<oo 

Let us introduce a function of y related to free energy 

fg )h (y) = -9[ln2 + ln®e,h(y)}- (20) 

To find the dependence of free energy on thermodynamic variables 9 and B one should find 
extreme points of fx,h (y)- F° r these points [dfK, y (y) /dy) 9 B = 0. Since 

jf] = - J K < \ [tanh (Ky + h) - y] <£> K , h (y) = -9 C [tanh (Ky + h) - y] , 
d VJe,B ®K,h{y) 

the state of minimum of free energy occurs for y obeying the equation 

y = tanh (Ky + h) . (21) 

For various values of 9 (K) and B (h) solutions of this equation provide the function y = 
y (9, B) of state variables. Therefore, the function $ (9, B) = max &K,h {y) is a composite 

— oo<y<oo 

function of 9 and B, namely §k, v (y {&, B)) and also depends implicitly on these two state 
variables. Now Eq. (|T9l can be rewritten in the form 

/ (0, B) = -9 ln 2 - 9 ln $ (9, B) , 

or 

/ (9, B) = -9\n2-9 ln j e -^ 2 (^)/2 cosh { Ky ^ fl) + fc] j . (22) 
Using two familiar identities^ dt&nhx/dx = cosh _2 x and 

cosh 2 x = (l - tanh 2 x) ^ = (l - y 2 ) _1 , (23) 
we calculate the second partial derivative of / with respect to y 



dy 



e.h 



= -9K{K[l-y 2 (9,B)]-l}. (24) 



In the Appendix we show that there exist solutions of Eq. (|2~T!) for which this derivative is 
positive, hence for them free energy is minimal. 

For the reversed magnetic field the solution of Eq. (1211) is — y 

-y = tanh [K (-y - h)} . (25) 

Therefore, free energy / (9, B) (I19p is an even function of B (and h) 

- HlZ^l = In 2 + In L-^-y^B)f/2 cogh { _ R { _ y ^ B)) _ h] \ = JJ^l . (26) 
9 I ) 9 



These properties of free energy and solutions of Eq. ([21]) are the reason why in the following 
we will always have in mind the positive value of B. 

The product of (dy/dh) B and (d 2 fK,h (y) /9y 2 ) dh is positive because 

(dy\ ( d 2 f K , h (y) \ 0K 

\dh) 9 \ dy* ) eh coBh 2 {Kv + h)' { ) 

Introduce here K c = 1 - the critical value of the parameter K. From the definition of K 
it is seen that the critical value of 9 is 9 C = J. 

In the mean field approximation (cf. for example^) the Hamiltonian ([3]) is approximated 

by 

N 
1=1 

For this Hamiltonian the statistical sum is 

Z mf (9,B) = [2cosh(Ky + h)] N 
and the corresponding free energy is 

F mf (9, B) = -9N [In 2 + In cosh (Ky + h)] . 
Using the thermodynamic identity (ITTj) we obtain Eq. (|2ip . 



IV. FREE ENERGY OF C-W MAGNET IN THE ABSENCE OF THE MAGNETIC 
FIELD 

Suppose that B = (h = 0). If we plot g (y) = y and q (y) = tanh (Ky) as functions of y, 
the points of intersection determine solutions of Eq. (12 ip . Referring to Fig. [1] (left panel), 
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we have to make distinction between cases. If 9 > 9 C (9 > J) the slope of of the function 
q (y) at the origin K = J/9 = 9 c /9 < 1 is smaller than the slope of linear function g (y) = y, 
which is 1, thus these graphs intersect only at the origin. It is easy to check that for this 
solution the second derivative of free energy f[T9"j) is positive (cf. Appendix). Therefore, the 
extreme indeed is a minimum (cf. Fig. [2]). 

On the other hand when 9 < 9 C (9 < J), the initial slope of tanh(Ky) is larger than that 
of linear function, but since values of tanh function cannot take values outside the interval 
(— 1,+1), the two functions have to intersect in two additional, symmetric nonzero points 
±y (6) (Fig. [1]). In this case in the Appendix we show that the second derivative of free 
energy is negative at the origin y — 0, which means that there is a maximum at y — (cf. 
Fig. [2]). This derivative is positive at y = ±y(9). Free energy / (8) attains minimal value 
at y = ±y (9), hence these solutions correspond to the thermodynamically stable states. 




FIG. 1. Graphical solution of Eq. (|2ip . The full line represents the function g(y) = y. The dotted 
lines: 6 < 9 C , the dashed lines: 6 > 9 C . Left panel: h = (B = 0). Right panel: h > {B > 0). 



The parameter K and temperature 9 can be expressed in terms of t — (9/9 c — 1), namely 



K 
9 = 



9c 
9 



1 



i + r 

i+t) 9 C . 



(28) 
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V. MAGNETIZATION AND MAGNETIC SUSCEPTIBILITY OF CURIE-WEISS 
MAGNET 



Consider the C-W magnet when the magnetic field B is brought back. Magnetization 
per particle, m, is a partial derivative of free energy / (9, B) after B 



m (9, h) 
d In § Kyh (y) 



\df (9,hY 




'0]n$(0,/i)~ 


dB 


= V 

e 


dh 



dh 



9 In <& K>h (y) 



dy 



y=y(e,B) 



dy(9,B) 
dh 



(29) 



y=y(o,B) 

Since $>K,h (y) attains an extremum at y (9, B), the second term on right hand side of (I2"9"j) 
vanishes. The contribution of the first term yields the equation of state 



K 

m — u tanh I — m + h 



(30) 



which has a closed analytic form. Therefore, according to Eq. ( 12"T|) . the study of y is 
equivalent to the study of of magnetization m. 

Consider the solutions of Eq. ( 12 lj) for B > 0. When 9 < 9 C the plots of functions 
9 (y) = V anc i Q (y) = tanh (Ky + h) intersect in three non symmetric and nonzero points 
(cf. right panel of Fig. [T]). Only for positive value of y = y(9,B) free energy / attains 
the global minimum. To one of negative values there corresponds a local minimum, to the 
remaining a maximum (cf. Fig. |2]). The negative values of y(9,B) (as well as m(9,B)) 
do not correspond to stable states and should be omitted. One should notice that negative 
values of y for positive h are not compatible with the symmetry (|25|) of the state equation. 
When 9 > 9 C the graphs of g (y) and Q (y) intersect at one point y (9, B) > 0. For this value 
of y (9, B) free energy has global minimum (cf. Fig. |2]). 

Consider Eq. (I2T|) for small values of \t\ and h = 0. We can expand tanh (Ky) into 
Taylor's series. Using Eq. fT28|) we obtain 

a/2 



y 



K 3/2 



1/2 



(31) 



Hence, in agreement with Fig. |3l we conclude that 



m ~ \t 



1/2 



(32) 



When positive B — > then y (9, h) — >■ y (9) of Fig |3j Thus, we note the existence of 
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e>e c , h=o 



e>e c , h=o.2 




FIG. 2. Dependence of function Jk,k (y) (|19 jl on y for paramagnetic (8 > 6 C ) and ferromagnetic 
(8 < 8 C ) phases in the magnetic field (h > 0) and when the magnetic field is absent (h = 0). 

spontaneous magnetization m (Q) = ny (Q) 

m (0,h = + ) = l ~ (33) 

yrriQ 6 C > 6. 

Even after turning off the magnetic field, below critical temperature the system remains 
magnetized, depending on the sign of B before its removal. The dependence of y for B = 
on 9 is singular. It is seen in Fig. [3] that at the point 8 = 8 C tangents of two branches of 
the curve are different. We note that the point 9 = 8 C is the boundary between the region 
of existence and nonexistence of magnetization, i.e. it is a critical point. We conclude that 
magnetization m is the order parameter. 

In the presence of even an arbitrarily weak (0 < /i < 1), the magnetic field magnetization 
m does not vanish below and above critical point (cf. Fig. [3] the right panel). The external 
magnetic field lowers the symmetry of the paramagnetic phase. From the point of view of 
magnetization the difference between paramagnetic and ferromagnetic phases vanishes and 
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the critical point ceases to exist. 





0.2 0.4 0.6 0.8 1 1.2 1.4 



0.2 0.4 0.6 0.8 1 1.2 1.4 



FIG. 3. Dependence of y on temperature. Left panel: h = 0. Right panel: h > 0. 
Consider the magnetic susceptibility 

xt (e,B) 



dm {9, B) 



OB 



(34) 



According to the definition ()34j) and Eq. (1291) the susceptibility is related to the second 
derivative of free energy with respect to induction 

'd 2 f(6,BY 



Xt (0,B) = - 



dB< 



(35) 



Differentiating both sides of Eq. fl30l) with respect to B and solving the obtained equation, 
we obtain the general expression for the susceptibility 

/i 2 [1 - y 2 (9, B)] 



Xt(9,B) q 1 _ K [ 1 _ y 2( 9)B )y 

Using Eq. fl24l the susceptibility can be written as 

l-y 2 (9,B) 



Xt {0,B)= l i i K 



fd 2 f K ,h(y)\ 




\ Q y 2 )o,h 


y=y(e,B) 



(36) 



(37) 



Since for a stable state the second derivative (d 2 fxh (y) /dy 2 )n h is nonnegative and 

y=y(9,B) 

< \y\ < 1, the susceptibility is nonnegative too. Notice that we succeed in linking together 
the macroscopic ( 1T3|) and microscopic ( 1361) stability conditions. 

For the paramagnetic phase {9 > 9 C (t > 0)) y — 0, and in the vicinity of the critical 
temperature we obtain 



Xt 



9 C (1+t) t 
12 



(38) 



This means that the critical index 7 = 1. 
The function arctanh?/ obeys the equation 



arctanhy = Ky + h . 



(39) 



We shall study solutions of Eq. (|39|) in vicinity of the critical temperature. For small y one 
can expand arctanhy into Taylor's series^. For 9 < 9 C (t < 0) Eq. fl39|) reduces to 



r-3 





t\ 




1 - 




\t\ 



y - 3h = 0. 



(40) 



In the ferromagnetic phase and in the absence of the magnetic field, the order parameter 
does not vanish 



3 1*1 



y 



(t < 0) . 



(41) 



When B = from Eq. fT4Tl) it follows that for ferromagnetic phase in vicinity of critical 
temperature m ~ a/3 \t\. Thus the critical index of magnetization is j3' = 1/2. 

Consider the susceptibility fl36|) in ferromagnetic phase in the vicinity of the critical 
temperature. Using the expressions ()28l) and (I4T1) for the ferromagnetic phase we obtain 

(42) 



and according to the definition §5§ for both phases the critical indexes of magnetization are 
equal: 7 = 7' = 1. We shall note that in the case of ferromagnetic phase one should use the 
relation (l4"0]) . We conclude that for both phases the magnetic susceptibility of Curie- Weiss 
magnet is divergent at the critical temperature. This singular behavior is shown in Fig. HI 



» o =l,h=0, [1=0.1 



e c =i,h=o.oi,n=o.i 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




0.2 0.4 0.6 0.8 1 1.2 1.4 




FIG. 4. Dependence of the susceptibility on temperature. Left panel: h = 0. Right panel h > 0. 
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At the critical point t — 0, hence y = y^ih 1 / 3 . Therefore, according the definition OH]), 
the critical index for the critical isotherm 5 = 3. 

Continuous phase transitions occur when a new state of reduced symmetry develops con- 
tinuously from the disordered (high temperature) phase. The ordered phase of M-C magnet 
has lower symmetry than the symmetry (j4j) of the Hamiltonian, thus the symmetry is spon- 
taneously broken. There exist two equivalent symmetry related states of C-W magnet with 
magnetization +m i — m respectively, with equal free energies. These states are macroscopi- 
cally different, so thermal fluctuations will not bring them into contact in the thermodynamic 
limit. To describe the ordered state we introduced magnetization - the macroscopic order 
parameter that describes the character and strength of the broken symmetry. 

VI. APPROXIMATE THEORY THE ANALYSIS OF ROOTS OF THE CUBIC 
EQUATION FOR MAGNETIZATION 

Now we shall study the roots of the cubic equation ( flOl) for the ferromagnetic phase. 
Since free energy is an even function of B we assume that h > 0. We shall study the so 
called the incomplete cubic equation 



where p = \t\ / (1 - \t\) , q = -3h/2. 

It is worthwhile to recall that assumptions proposed by Landau for an incompressible 
magnet result in free energy depending only on even powers of magnetization^^ 1 ^ 



y 3 + 3(-p)y + 2q = 0, 



(43) 



/ (m, T) = / (T) + a (T) m 2 + ^ (T) m 4 . 



This form of free energy leads to a cubic equation for magnetization. 
Introduce a characteristic value of the parameter ht 




(44) 



Roots of Eq. ()43l) depend on the sign of the discriminant 15 D = (q 2 





(45) 



If D < inequalities —h t < h < h t hold. If D > 0, then h > h t or h < 



h t . 
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If D < all tree roots are real 



y [ <] (t, h) = u< a, h) + v< (t, h) , 

yi <] (t, h) = e 2U< (t, h) + e lV< (t, h) , (46) 



where 



(47) 



u< {t, h) = -q + iy/\E\, w< (t, h) = [u < (t, h))* , 
e x = + /2, e 2 = e\. 
Using relations (j4"7]) we can show that 

y [< (t,h) = 2Rew< (t,h), 
(yl <) Y(t,h)=y^(t,h) (a = 2, 3). 

If D > (i.e. /i 2 > /i 2 ) only one root is real 

y> (t, h) = U> (t, h) + V> (t, h) , (48) 

where 

u> (t, h) = \Jy + vb > v > (t ' h) = \Jy~^- (49) 

The two remaining two roots are complex. 

If we combine Eqs. (J4"4l - (j4"9|) we obtain the functions u K (t,u), v < (t,u) and u> (t,u), 
v> (t, u) in the useful form 



(t,h) = ( 3 /2) 1/3 ^77^^f, v > (t,h) = (3/2) 1/3 ^- v/^^f, 
u< (t, h) = (3/2) 1/3 \f h + i^h 2 -h 2 , v K (t, h) = (3/2) 1/3 ^fh-i^/h 2 -h 2 . 

Consider functions u > (t,h) and v > (t, h) for negative h 



(50) 



u> (t, -h) = ^=1 \Jh- \/h 2 - h 2 = V^Tvy (t, h) , 

I ( 51 ) 

v> (t, -h) = S^I y h + \fh 2 ^h 2 = ^1 U> (t, h) . 

Similar relations hold also for u K and v < . Among three root of unity £i, £2 and £3 = — 1 
only the latter root yields proper symmetry relation (|25|) 

(f , -fc) = -yf (t, h) (a = >, <; j = 1, 2, 3) . (52) 
15 



In Fig. [5] we have plotted the dependence of roots (|46|) and (jUJ) on the magnetic field 
h. For positive values of and —ht < h < 0, as well as for negative values of y< 
and < h < h t , the signs of h and these two roots do not agree. The symmetry (1251) of 
equation of state is broken. On the line BB' (corresponding to the root y< ) the derivative 
(dy (6>, h) /dh) e < is negative. From Eq. (l2~Tj) it follows that (d 2 fe^ h (y) /dy 2 ) e h is negative 
too. In this interval of values of h free energy has maxima, therefore CW magnet is not in 
stable states. We conclude that for the root y< ' (i, h) free energy attains maximal values. 

" I 7" 1 n ! 1 — 1 1 



-0.4 




-0.6 -0.4 -0.2 0.2 0.4 0.6 



FIG. 5. Plot of dependence of roots of Eq. (T33|) on h for h t = 0.1. For \h\ > h t the curve represents 
the function y> (t, h). In the interval —h t < h < h t the plot consists of (t, h) (m = 1, 2, 3) 



(3) 

y< 



For < h < h t the root ?/< is positive. The remaining two roots are negative and 
< y< (cf. Fig. [6]). Further y^ — > when /i — >■ = h t — e; (0 < e « 1) and 



( dy®(e,h) , 
lim I = — oo. 
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FIG. 6. Dependence of 
From Eq. fl27]) it follows that 



(3) 

y< 



and y< on h. 



lim 

h—¥ht~ 



d 2 fe )h (y) 



0. 



9y 2 / e,h 

This means that for h = ±h t one should consider the third derivative of free energy with 
respect to y. 

(3) 

In the interval —h t < h < h t the derivative of the root y K is positive, therefore the 
derivative (d 2 f e<h (y) /dy 2 ) 



9,h 



(3) 

v=y< 



is also positive. This means that in this interval free 



energy has a (local) minimum. Free energy exhibits the global minimum for the root y< 
for h t > h > 0. This behavior of free energy and of solutions of cubic equation is shown in 
Fig. ED 

For h > ht and h < —ht free energy has the global maximum for the real root y > (t, h) 
(cf. Fig. ED. 

Note that as a result of approximation yielding Eq. (140 p values of \y\ may exceed the 
limiting value 1. This means that we shall restrict ourselves to small values of \t\ and h t . 

We shall point out that value h t (|44p of parametr h define such value B t of magnetic 
induction B, for which the value of induced magnetization m in d ~ xt • B t (B t = 9 c h t / ft with 
h t given by Eq. fl38l) ). is in accordance with Eq. fj4~Tl) . 
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If h <C h t (\t\ 7^ 0) the magnetic field B is week and does not influence the thermodynamic 
quantities characterizing the system. If h ^> h t the field B is strong. If t — (T — T c ) all 
magnetic fields are strong. As we have shown, if t ~ and field is strong, m ~ 



VII. PROPERTIES OF THE INTERNAL ENERGY, ENTROPY AND SPECIFIC 
HEAT OF THE CURIE- WEISS MAGNET 

To find the internal energy U we shall use the familiar thermodynamic identity 11 

For the internal energy per one spin this formula gives 

u = -e 2 ^- \-Ky 2 /2 + In cosh (Ky + h)] . 

As a result of simple calculations we obtain 

u{6,B) = - J -y 2 {6,B)-^By. (53) 

The first term of this equation is the interaction energy per one spin, whereas the second 
term is the energy of a spin in the magnetic field. When 9 = 0, y = 1, hence u = —J/2 — Bfi. 
When B = 0, in the paramagnetic phase y = and u = 0. When B = in the ferromagnetic 
phase in vicinity of 6 C one has K ~ (1 — and 



y « ±v/3|t|/ (l-|t|), (54) 

thus, m ~ |i| . As we can see from Fig. [3, in the absence of the magnetic field the behavior 
of internal energy is singular at K = 1. The curve representing the function u(9,h = 0) 
consist of two branches. At 9 = 9 C their derivatives are different. When the magnetic field 
is turned in this singularity is washed out. 
Consider entropy 

Calculating the derivative we obtain for entropy per particle the familiar thermodynamic 
identity^ 

s(9,B) = ^-[-f(9,B) + u(9,B)]. (56) 
18 
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FIG. 7. Dependence of the internal energy on temperature. Left panel: h = 0. Right panel h > 0. 
For low temperatures 



f(0,B) 



u{9,B) 
9 



Therefore, even when the magnetic is turned in, at 9 = entropy vanishes. The spins are 
completely ordered, therefore entropy acquires the lowest value. Since in the absence of the 
magnetic field the internal energy is singular at K c , entropy is also singular at the critical 
point. When B = in the paramagnetic phase (9 > 9 C ) y = 0, and the internal energy 
vanishes. From the definition ( !T9|) it follows that —ksf (9 > 9 C , B = 0) /9 = /c£ln2, and 
s = ksln2 ~ 0.7 x ks- Spins in the paramagnetic phase are completely disordered and 
entropy reaches its greatest value. The dependence of entropy on temperature is shown in 
Fig. E 
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FIG. 8. Dependence of entropy on temperature. Left panel: h = 0. Right panel: h > 0. 

If B = in the paramagnetic phase y = 0, hence u (8 > 8 C , h = 0) = 0. In the fer- 
romagnetic phase and in vicinity of 9 C according (13TT) u ~ \t\ . For 9 = 9 C the internal 
energy vanishes u = 0. This means that the internal energy is the continuous function of 
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temperature u(9~) = u(9+) = 0. 

Entropy s also is a continuous function of temperature. To show this property in the case 
of ferromagnetic phase and vicinity of 9 C we use Eqs. (J22J) and ( 1ST]) and Tylor's series for 
In [tanh (Ky)] we get 

s « A; B (In 2 - 3 |t| /2) . 

We see that, as one may expect, that entropy of the ferromagnetic phase is smaller than 
entropy of paramagnetic phase and for \t\ = attains its maximal value. 

The behavior of heat capacity in vicinity of the critical temperature is more complex. To 
calculate heat capacity per one particle one can use one of two thermodynamic relations, 
namely^ 

c(9 , BHfe wn , (57) 



89 

or 

'8s (9,B 



B 



'(^"'i-irj, (58) 

It is an easy task to show that these identities yield the same result. We shall consider Eq. 
(!57p . From Eq ( |5B1 it follows that heat capacity c(9,B) depends on derivative (8y/89) B 

c = -k B e(Ky + k) 

Differentiating both sides of Eq. ( 12 ip with respect to y, solving the obtained equation for 
(8y(9,B) /d9) B and applying the identity (l2~3"j) . we obtain an analytic expression for this 
derivative 

'8y\ 1 {Ky + h) 



89 J B 9 cosh 2 (Ky + h) — K' 
With the help of the above relation we obtain the final form of expression for heat capacity 
per one particle 

(1 - y*) (Ky + hf 
C=kB (1-K) + Ky* ■ (59) 

When B = in the paramagnetic phase (9 > 9 C ) magnetization m vanishes, i.e. y = 0. 
Hence, c = 0. Since limy = 1 in this limit c also vanishes. 

Consider heat capacity in ferromagnetic phase in the vicinity of the critical temperature. 
Using expressions ( |28i) and ( 15^1) we obtain 

£»§<i-H)- («» 
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Heat capacity is discontinuous at 6 C 



c (0- B = 0) - c (0+ ,6 = 0)= 3A; B /2, (61) 

with 0± = C ± e, and < e < 1. The zero field heat capacity is singular at the critical 
temperature, whereas in the magnetic field it exhibits a peak at the transition point (cf. 
Fig. E]). 



e c =i, h=o e c =i, h=o.oi 
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FIG. 9. Dependence of heat capacity per one particle on temperature. Left panel: h = 0. Right 
panel: h > 0. 



VIII. NUMERICAL RESULTS 

Until now we described the phase transition using the approximate expression for 
arctanh (Ky + h). Now we shall give up this approximation, which means that we shall rely 
on numerical calculations. Such approach will shed light on the approximation used in Sect. 

IYU 

Let us rewrite Eq. in the form 

arctanhy — - - = h. (62) 

For a given value of \t\ we numerically solve Eq. (1621) for various values of h. We plot 
y = yt (h) in Fig. [TUJ The obtained plot resembles the plot obtained for roots of Eq. (|43|) . 
However, the critical values of the parameter fi^ num ^ are smaller than obtained h[ approx ^ as 
a result of calculation based on the approximate theory presented in the Sect. |VH Besides, 
unlike for the approximate theory, values of y num are confined in the interval [-1,1]. In Fig. 
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FIG. 10. Dependence of y% (h) on the magnetic field h. Dashed line represents the dependence of 
roots of Eq. (|43p on the magnetic field. The full line is plot of the function (|62p . 

□J] we compare values of h { t approx) and hf um) . It is seen that for small values of |t| both plots 
differ a little. With growing \t\ difference is more pronounced. 



IX. APPENDIX 



We shall show that there exist solutions y (9, B) of Eq. (|39|) for which the second derivative 
of function f K>h (y) (Eq. f l2"Uj) ) with respect to y calculated at y — y (9, B) 



d 2 f 9 , h (y) 

dy 2 



e.ii 



= -9K{K[l-y 2 (9,B)] - 1} (63) 

y=y(8,B) 



is positive. 

Consider temperatures higher than 9 C . 

1. For B = in the paramagnetic phase K < 1 and magnetization vanishes. In this case 
y = 0, therefore, the second derivative f )63|) [— 9K (K — 1)] > is positive. 
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FIG. 11. Plot of dependence of critical value of the magnetic field ht on the parameter t. The full 
line - result numerical calculation, the dashed line represents the dependence resulting from the 
cubic equation (|43j) . 



2. If B > magnetization is positive, hence 1 > y > 0. The double inequality < 
K [1 - y 2 (9, B)} < 1 holds. Therefore, K [1 - y 2 (6, B)] - 1 < 0. Hence, the the 
second derivative (1631) is positive. The same arguments are valid for the negative 
value of induction (B < and y < 0). 



In the case of 6 < 9 C the parameter K is greater than unity. As we know (cf. Sect. [13 
if B = 0, Eq. fT2Tl) has three solutions, namely y = and y = ±y (6) = ±2/ . 



1. For y = the derivative (l63l) is negative 

'd 2 fe, h (y) 



dy 2 



e.ii 



6K{1-K) < 0. 



(64) 



This means that in the ferromagnetic phase the solution y = corresponds to a max- 
imum of free energy. 



2. In the case of two remaining solutions the second derivative reads 

'd 2 f K , h ( y y 



dy 2 



e,h 



-9K [K (1 - yl) - 1] 



(65) 



y=±yo 
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We shall express the parameter K by yo- Using Eq. (|39|) we can write 

arctanhyo 



K -- 

With the help of this relation we find 

'd 2 f K , h ( y y 



yo 



dy 2 



e,h 



-9K 



y=±yo 



arctanhyo 
Vo 



(l-^/o 2 ) 



(66) 



For small y we can use the approximate expression arctanh?/ ~ yo + yo/3- This yields 
the inequality 

(1 — ?/q) arctanh?/ 



(l +2/0 2 /3) (1 - yl) < (1 - y 4 ) < 1, 



and the second derivative of free energy is positive. In the ferromagnetic phase, in 
vicinity of critical temperature, for solutions ±?/o free energy reaches a minimum. 

The function 

arctanhy , 2 x 

<f{yo) = I 1 - yo) 

yo 

is monotonically decreasing in the interval (0, 1] and < ip (yo) < 1. 

Using the logarithmic representation^ arctanhyo = \ In fjj~i) an d the limiting 
value 18 lim x In x = we can show that 

x^O 

arctanlw/o h 2 \ n 
km (l-y ) =0, 

vo-*i- Vo 

hence, for low temperatures the right hand side of Eq. fl66|) is positive. The plot 
of function (1661) is shown in Fig. [12] We conclude that in the ferromagnetic phase 
minimums of free energy occur for y = ±yo- 

3. Assume that B > (h > 0) and 9 < 9 C (K > 1). In this case, according to Sect. |V] 
Eq. ([21]) has three solutions yf ] (t, h) > 0, y ( 2 <] (t, h) < 0, y[ <] (t, h) < 0. Proceeding 
as before we express K by y[ K ^ and h 

arctanhy| < ' > — h 



K 



Now, Eq. ( |63j) takes the form 

'd 2 f K , h (yy 



dy 2 



-9K 



y[ <] 



arctanhyj^ — h 



(67) 



(<) 



Vi 



(<) 



i-ly[ <] 



n. 



(68) 
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df K , h (y)/dy 



d2f K,h(y) /d y 2 
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FIG. 12. Plot of /fc^ (y) (solid line), dfx,h (y) /dy (broken line) and d 2 fx,h (y) /dy 2 (dotted line). 
To zeros of dfx,h(y) /dy there correspond extremums of fx,h{y)- Two of them are minimums 
because d ^ K d £ {v) > 0. 



Since < < 1, the inequality 
of Eq. ( )68|) is non-positive 



i - (y{ <] 



h 



(<) 



i - (^ 



< 1 holds. Thus, the additional term 



<0, 



and, even in this case, the second derivative of free energy is positive. 
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